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In this study, considering the long-range interaction with an inverse-square and its trigonometric 
and hyperbolic variants in SCM model we investigate entanglement in (1/2,1) mixed-spin XY model. 
We also discuss the temperature and magnetic field dependence of the thermal entanglement in this 
system for different types of interaction. The numerical results show that, in the presence of the 
long-range interactions, thermal entanglement between spins has a rich behavior dependent upon 
the interaction strength, temperature and magnetic field. Indeed we find that for less than a critical 
distance there are entanglement plateaus dependent upon the distance between the spins, whereas 
above the critical distance the entanglement can exhibit sudden death. 
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I. INTRODUCTION 

Over the last decade or so, quantum entanglement has 
been recognized as crucial in various fields of quantum in- 
formation [l[ such as in quantum computation Q , quan- 
tum teleportation Q, superdense coding Q, quantum 
communication [5j, cjuantum perfect state transfer [6|, 
quantum cryptology 0, H| and quantum computational 
speed-ups [3, [l(| • Potential applications of entanglement 
in these fields have stimulated research on methods to 
quantify and control it. In the solid state arena, consid- 
erable attention has been devoted to interacting Heisen- 
berg spin systems. Spins are recognized as candidates 
for solid state quantum computation [lll - fTBj and even 
short range communication [lq . \vf\ . Interactions such as 
Heisenberg enable gate operations in solid state quantum 
computation processors [ldl - tl5j . Therefore significant re- 
search has been performed to understand quantum en- 
tanglement behavior in spin systems, such as all various 
kinds of Heisenberg (XX, XY, XXZ and XYZ) models, 
and similarly Ising models. Many detailed and extensive 
investigations have addressed the dependence of entan- 
glement on parameters such as magnetic field, single-ion 
anisotropy, temperature and interaction strength [l8l - l42] . 

Recently, a new typ e of long-range interaction has been 
used in Refs. |43l445| to obtain long-distance entangle- 
ment in the spin system. In these works, spin pair in- 
teraction is given with a factor inversely proportional to 
a power of the distance between sites as J (R) ~ R~ a . It 
is shown in these studies that long-distance entanglement 
can be obtained using this interaction type for different 
values of the interaction parameter a in the Heisenberg 
spin systems. Here we must remark that this interac- 
tion type and its variants except the entanglement have 
been considered to explore many physical phenomena. 
Indeed, the inverse-square, trigonometric and hyperbolic 



interacti ng p article systems [46l - |48l ] and its spin general- 
izations (49l - l52| are important model of many-body sys- 
tems due to their exactly solvability and intimate connec- 
tion to spin systems in condensed matter |49|-|59| , and in 
the other areas in physics. These interaction types are 
called different names such as Haldane-Shastry models, 
Calogero-Moser or Sutherland-Calogero-Moser models or 
interactions in the literature. But, it is assumed that gen- 
eralized Haldane-Shastry models is the supersymmetric 
partners of the Sutherland-Calogero-Moser type models. 
Hence, for historical reasons, in this study we will express 
these interactions as Sutherland-Calogero-Moser (SCM) 
model or SCM tpye interactions. 

In this study, inspired by Refs. [H|-|45j], we will con- 
sider long-range interaction with an inverse-square and 
its trigonometric and hyperbolic variants given in the 
SCM model to obtain entanglement between two spin in 
the mixed-spin (1/2,1) XY Heisenberg spin system. Our 
numerical results show that the long-range interactions 
lead to topological plateaus in entanglement. That is, 
the entanglement between spins still subsists throughout 
plateau up to critical distance although spins are sepa- 
rated from each other. 

This paper is organized as follows. In Section II we de- 
fine the mixed spin-1/2 and spin-1 XY model, and give 
the eigenvalues, eigenvectors, density matrix of the sys- 
tem and definition of negativity as a measure of entan- 
glement. In Section III, we present and discuss the nu- 
merical results of modeling, giving the negativity of the 
system for SCM types i.e. inverse-square and its trigono- 
metric and hyperbolic variants with varying model pa- 
rameters. Finally Section IV is devoted to conclusions. 



II. MODEL 
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The Hamiltonian for the isotropic mixed spin (1/2,1) 
Heisenberg XY spin chain under an applied magnetic 
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field B is given by 



H = ^B(st 
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where J(R) is the spin interaction coupling which will 
be defined in terms of SCM interactions. Also Si and Si 
(i = 1,2) are spin operators of the spin-1/2 and spin-1 
components, respectively. The basis vectors of the two- 
spin system defined in Eq. (1) are { | — 1/2,— 1), 1 1/2, 0} , 
|-1/2,1), |l/2,-l), |-l/2,0), |1/2,1) }, where \a,S) is 
the eigenstate of s z and S z with corresponding eigenval- 
ues given by s and S, respectively. The corresponding 
eigenvalues and eigenvectors of H for all J(R) can be 
given as follows 



H\ 



-1/2,-1) = -3/2B|-l/2,-l) 
F|1/2,1) = 3/2B| 1/2,1) 
H\^) = W± 1^) 



H |V' ± ) 



where 



^ = 1/^2(11/2,0) ±|l/2,-l)) 
^± = 1/^2(1-1/2,1) ±|-l/2,0)). 



(2a) 
(2b) 
(2c) 
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Here (f> and ij) are maximally entangled states — so- 
called Bell states. The eigenstates of the isotropic XY 
model given by Eq. (2) with Eq. (3) do not change for 
different J{R) values. However, the eigenvalues W± and 
Q± in Eq. (2) do change dependent upon the value of 
J(R). Each eigenvalue, for different types of the SCM 
interaction, will be explored in detail in the next section. 

The partial transpose of the density matrix p for the 
two-spin which is defined by Eq. (1) is written in the form 



/an 012 \ 

0,21 a 22 

a 33 

a 44 

a 55 a 56 

V a 65 a 66 / 



(4) 



where a,-j = 1, 6) are matrix elements of p Tl . For 
different cases of J, the matrix elements of p Tl are pre- 
sented, along with the negativity. 

Negativity is a measure of the quantum entangle- 
ment appropriate for application to higher spins [60H62j . 
Therefore we will use the concept of negativity to study 
entanglement in (1/2,1) mixed-spin XY systems. The 
negativity can be obtained using the density operator 
of the quantum system. The state of a system at ther- 
mal equilibrium can be described by the density operator 
p (T) = exp (-/3H) /Z, where Z = Tr (exp (-f3H)) is the 
partition function and /3 = 1/kT (k is Boltzmann's con- 
stant, which is set to unity k = 1 hereafter for the sake 
of simplicity, and T is the temperature). By choosing 
an appropriate set of orthonormal product basis states 



{\viVj)} = {\vi) \vj)} for the density operator, the partial 
transpose is defined by its matrix elements 



Pmn,nv ~ ( v mVfj, \ p \v n Vv) — Pmv;i 



The negativity of a state p is by definition 



(5) 



(6) 



where \ii is the negative eigenvalue of the partial trans- 
pose density matrix p Tl . Here p Tl is the partial transpose 
with respect to the first system. The negativity N is re- 
lated to trace norm of p Tl via 



N 



(7) 



where the trace norm of p Tl is equal to the sum of the 
absolute values of the eigenvalues of p Tl . 



III. RESULTS AND DISCUSSIONS 

In this section in order to discuss entanglement we will 
compute the negativity to quantify entanglement for two- 
spin in the isotropic mixed-spin (1/2,1) XY Heisenberg 
spin chain given by Eq. (1) considering different long- 
range interaction types given in the SCM model [46l — 
[52[. Therefore, in the below Subsections, we will con- 
sider inverse-square interaction coupling Jq/R 2 as type 
I, trigonometric interaction coupling Jo / sin 2 (R) as type 
II and hyperbolic interaction coupling Jo/sinh 2 (i?) as 
type III in Subsections III.A,B and C respectively. 



A. SCM Type I: J(R) = J /R 2 

Firstly we consider the inverse-square interaction type 
in SCM model, which is defined with exchange interac- 



tion J(R) = J /R 2 [46H52J- Here we set J = 1 for sim- 
plicity. The eigenvalues for this type of SCM interaction 
are given by 



Using the eigenvalues in Eq. (8) and the eigenvectors in 
Eq. (2) (with Eq. (3)) the density matrix p can be con- 
structed. For a SCM interaction of type J(R) = l/R 2 , 
the partial transpose of p shown in Eq. (4) has matrix 
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elements given in Eq. (9) as follows: 



a 12 = a 21 =--e- sinh[^^] 
a 22 = a 33 = ie-*co S h[^^] 
a 44 - a 55 - cosh[ ^ R2J1 ] 

a 56 = a 65 = -i e -*smh[^=^;] 
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where the partition function Z is written as 



Z = 2{ cosh[^] 



B 1 
2 coshf— — 1 cosh[^= 



(10) 



For the SCM interaction of type J{R) = 1/R 2 , the nega- 
tivity given by Eq. (7) for the two-spin system defined in 
the Heisenberg mixed XY model in Eq. (1) can obtained 
from summation of the negative eigenvalues of the partial 
transpose matrix Eq. (4) with elements as in Eq. (9). 

The negativity of the two-spin system in the case 
J(R) = 1/R 2 is plotted versus R for different temper- 
ature T values at fixed magnetic field B = 1 and for dif- 
ferent magnetic field values at fixed temperature T = 1, 
in Fig. (la) and Fig. (lb), respectively. Firstly we note 
that there is a singularity at R = 0, hence, as it can 
be seen from both figures that negativity drops to zero 
in the limit R — > 0. On the other hand, as it can be 
seen from Fig. (la), the negativity reaches a plateau for 
small R at a fixed temperature T and fixed magnetic 
field B. The value of the negativity plateau is 0.5 at very 
low temperatures, decreasing at higher temperatures al- 
though with the plateau behavior still apparent. Similar 
behavior is observed in Fig. (lb), with the negativity at- 
taining a plateau at small R. For relatively high magnetic 
fields the plateau value approaches 0.5, but decreases for 
lower magnetic values as expected. A very interesting 
aspect of the behavior in these figures is that the neg- 
ativity shows this plateau behavior, like in topological 
quantization, dependent upon the distance between the 
two spins for different temperatures at fixed B and sim- 
ilarly for different magnetic fields at fixed temperature. 
It is interesting to note that the plateau width (in R) 
increases as it approaches its maximum for lowering T, 
whereas the width decreases as it approaches its maxi- 
mum for increasing B. In both figures, for increasing i?, 
in some cases the negativity exhibits a form of sudden 
death (63l - l65| . reaching zero at a critical R value of R c 
[6r| . The results of Fig. (la) demonstrate the dependence 
of the critical R c on temperature T at fixed B, whereas 
the results of Fig. (lb) demonstrate some insensitivity of 
R c to magnetic field B at fixed T but with the entangle- 
ment slope (at i? c , as a function of R) dependent upon B. 



FIG. 1. Negativity as a function of R for the two-spin fer- 
romagnetic XY system (a) for selected temperatures at fixed 
magnetic field B = 1, (b) for selected external uniform mag- 
netic fields at fixed temperature T = 1. 
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FIG. 2. Magnetic field dependence of negativity in the two- 
spin ferromagnetic XY system for selected temperatures (a) 
at fixed R = 0.5, (b) at fixed R = 1.5. 



Here we express that the topological plateaus in Fig. (1) 
emerge from the presence of the long-range interaction 
with inverse-square 1/R 2 in this model. 

In order to further investigate the magnetic field de- 
pendence of the entanglement for different R values in 
the larger R regime, the negativity is plotted as a func- 
tion of magnetic field for different temperature values 
(T = 0.05,0.15,0.25), at fixed R = 0.5 in Fig. (2a), and 
at fixed R = 1.5 in Fig (2b). It can be seen from these 
figures that although the negativity curves have differ- 
ent detailed characteristics, the negativity for the two 
different values of R tends to zero for increasing B val- 
ues. Identifying a critical B c beyond which there is no 
entanglement suggests that this critical B c value and the 
manner in which the entanglement vanishes both depend 
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FIG. 3. Temperature dependence of the negativity in the 
two-spin ferromagnetic XY system for selected magnetic field 
values (a) at fixed R = 0.5, (b) at fixed R = 1.5. 

upon the chosen R and the temperature T. 

Similarly, to further investigate the temperature de- 
pendence of the entanglement for different R values in the 
larger R regime, the negativity is plotted versus temper- 
ature for different magnetic field values (B — 0.5,2,4), at 
fixed R = 0.5 in Fig. (3a), and for different magnetic field 
values (B = 0.05, 0.15, 0.35) at fixed R = 1.5 in Fig. (3b). 
Since the negativity does not appear for B = 0.5,2,4 at 
R = 1.5, we plot Fig. (3b) for different B values. As it can 
be seen from these figures, the temperature dependence 
of the negativity for different R and B values has very 
similar characteristic behavior. In both figures, the neg- 
ativity curves for different magnetic values meet at the 
same temperature point which indicates a critical tem- 
perature value T c . Furthermore, Fig. (3) demonstrates 
critical temperature T c dependence on R, consistent with 
the critical R c dependence on T of Fig. (la). Similarly to 
Fig. (1) singularity behavior in negativity appears in the 
limit R -> in Fig. (3). 

To further demonstrate and unify the dependence of 
the entanglement on temperature and magnetic field, a 
surface plot of the negativity is given in Fig. (4), for an 
example fixed R = 0.5, as a function of both T and B. 
As can be seen from Fig. (4) , the negativity of the two- 
spin XY system has a plateau for this fixed R value. 
We also know from the previous Figs. (1), (2) and (3) 
that the range of the plateau depends on magnetic field 
and temperature and we can qualitatively comment that 
range (in T and B) of the plateau for small R values 
will be bigger than that for large R values. Furthermore, 
Fig. (4) shows clearly that for very low temperature val- 
ues the negativity exhibits critical behavior (coming off 
the plateau) and decreases rapidly with increasing field 
at a critical magnetic field value. However, for higher 
temperatures beyond T c when there is no access to the 
plateau, the negativity simply decreases smoothly for in- 
creasing magnetic field. 




1.5 4 

FIG. 4. Surface plot of the negativity in the two-spin ferro- 
magnetic XY system as a function of temperature and mag- 
netic field at fixed R — 0.5 



To summarize, in this subsection, we have investi- 
gated the behavior of the entanglement of a two-spin 
XY system in the case of long-range SCM interaction 
J(R) = 1/R 2 . We have discussed the dependence of the 
negativity on interaction parameter R, magnetic field B 
and temperature T. We have shown that entanglement 
appears in the presence of the long-range interaction in 
some parameter regimes, and there is critical behavior in 
entanglement which suddenly deaths at critical values of 
interaction parameter, magnetic field and temperature. 



B. SCM Type II: J(R) = Jo/sin 2 R 

Secondly we consider trigonometric version of long- 
range interaction in SCM model, which is defined with 
exchange interaction J(R) = Jo/ sin 2 R \Idt^5z\ . Here we 
set Jo = 1 and 8 = R for simplicity. The eigenvalues for 
this type of SCM interaction are given by 



±3B/2 , W± . ±^±V2, Q± _ TSazM. 

zip zip 



where ip = 3 — 4 cos[2i?] + cos[4R]. Using the eigenvalues 
in Eq. (11) and the eigenvectors in Eq. (2) (with Eq. (3)) 
the density matrix p can be constructed. For the SCM 
interaction type of J(R) = 1/ sin 2 R, the partial trans- 
pose of p shown in Eq. (4) has matrix elements given in 
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Eq. (12) as follows: 

an = -^e*r (12a) 

ai2 = 02i = f^(-e" 51 + e- 52 ) (12b) 

022 = o 33 = g(e 5l +e fe ) (12c) 

o 44 = a 55 = + e- &2 ) (12d) 

a56 = o 6 5 = g(e 51 -e 52 ) (12e) 

a 66 = ~e~%T (12f) 

where = ^os^csc^ 5 2 = and 

_ , csc[_R] 4 (3S+4Bcos[2_R]+-Bcos[4fi]+8\/2^/sm[fl] 4 ) . 

MT = ex P("F i6T /• 

The partition function Z for the interaction type J = 
l/sin 2 (i?) can therefore be written as 

Z = ,i_e-( 5l +^ V 1+2 ' 2 + e 2Sl+S2 + M+e^ (13) 
+M+e fe + {i+e^ 1 ^ + fi 2 + e Sl } . 

For the SCM interaction of type J(R) = l/sin 2 i?, the 
negativity given by Eq. (7) for the two-spin system de- 
fined in the Heisenberg mixed XY model in Eq. (1) can 
obtained from summation of the negative eigenvalues of 
the partial transpose matrix Eq. (4) with elements as in 
Eq.(12). 

The negativity for this J(R) = 1/ sin 2 R case is plot- 
ted versus R (0 < R < ir) for different temperature T 
values at fixed magnetic field value B in Fig. (5a), and 
for different magnetic field values at fixed temperature T 
in Fig. (5b). One can see that negativity drops to zero in 
the limit R — > and R — > tt in both figures. In these lim- 
its, the negativity have two singularities. On the other 
hand, it can be seen from these figures that the nega- 
tivity drops suddenly to zero with R, creating a valley 
symmetric about R — tt/2. This occurs for a fixed B = 1 
value at the various temperature values used for Fig. (5a), 
with the valley width and side plateau heights dependent 
upon the chosen T value and lower T creating a narrower, 
higher-sided (N = 0.5) valley. The effect also occurs for 
a fixed T = 1 value at the various magnetic field values 
used for Fig. (5b), with the valley width and side plateau 
heights dependent upon the chosen B value and higher B 
creating a slightly wider, steeper-sided and higher-sided 
(approaching N — 0.5) valley. Given the symmetric be- 
havior of the negativity in the interval < R < tt, as 
shown in Figs. (5) it is clear that there are two critical 
R c values in these figures, equidistant from R = tt/2. 

In order to further demonstrate and unify the depen- 
dence of the entanglement on R (0 < R < tt) and T, the 
negativity surface is plotted as a function of R and T for 
a fixed value B = 1 in Fig. (6). This surface plot shows 
the temperature dependence of the interesting valley be- 
havior of the negativity, with side plateaus, as a function 
of R. It is clear that at this fixed value of B = 1, the neg- 
ativity valley (symmetric in R about R = tt/2) widens 
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FIG. 5. Negativity as a function of R for the two-spin fer- 
romagnetic XY system in the case J(R) = 1/ sin 2 R (a) for 
selected temperatures at fixed magnetic field B — 1, (b) for 
selected external uniform magnetic fields at fixed temperature 
T — 1. 
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FIG. 6. Surface plot of the negativity for the two-spin ferro- 
magnetic XY system in the case J(R) = 1/sin 2 R as a func- 
tion of temperature T and interaction parameter R at fixed 
5 = 1. 

with increasing temperature, with the side plateaus di- 
minishing and dropping from their low T maximum of 
N = 0.5. The T-dependence of the two symmetric criti- 
cal R c values is clearly visible in this surface plot. How- 
ever, whilst the valley profile of the negativity with R per- 
sists all the way down in T to zero temperature, it is clear 
that at the lowest temperatures the valley floor exhibits 
a ridge for which the negativity is non-zero for all val- 
ues of R including the mid-point R = tt/2. At these low 
temperatures critical behavior with N dropping suddenly 
to zero thus disappears, although the negativity still ex- 
hibits a critical change from the plateau at N = 0.5 to 
the small finite value at the valley floor. The behavior of 
this low-T ridge in negativity as the temperature actually 
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FIG. 7. Surface plot of the negativity for the two-spin ferro- 
magnetic XY system in the case J(R) = l/sin 2 i? as a func- 
tion of magnetic field B and interaction parameter R at fixed 
T = 0.001. 
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critical R c values is clearly visible in this surface plot and 
it is also clear that these merge at a specific B point and 
below this value of B the entanglement is at its maximum 
for all values of R, with no critical behavior. 

The final figure for the SCM interaction J(R) = 
1/ sin 2 R is Fig. (8), which shows the behavior with B and 
T of the negativity surface at a fixed value of R = tt/2, 
which is the mid-point of the R range — about which the 
valley behavior in R in Fig. (6) and Fig. (7) is symmetric. 
It can be seen that non-zero negativity appears at low T 
for this fixed R — tt/2, consistent with the ridge seen in 
Fig. (6). The behavior of the ridge (or valley floor) mid- 
point as a function of B and T is captured in Fig. (8). It 
can be seen from this that at very low and high magnetic 
field values, the negativity (and thus the ridge) disap- 
pears as T —¥ or may not even appear at all, whereas 
for intermediate B ~ 1/2 the negativity (ridge) appears 
and grows as T — ► 0. 

As a result, in this subsection, we have investigated 
the behavior of the entanglement of a two-spin XY 
system in the case of the long-range SCM interaction 
J(R) = l/sin 2 i?. We have discussed the dependence of 
the negativity on interaction parameter R, magnetic field 
B and temperature T. We have demonstrated the sym- 
metric behavior about R = tt/2 (with singular behavior 
as R —> and R — > tt). The critical plateau and valley 
behavior is seen in significant parameter regimes along 
with the associated entanglement sudden death. 



C. SCM Type III: J{R) = J /sinh 2 R 



FIG. 8. Surface plot of the negativity for the two-spin ferro- 
magnetic XY system in the case J(R) = l/sin 2 i? as a func- 
tion of temperature T and magnetic field B at fixed R = tt/2. 



approaches zero is dependent upon the fixed value of B 
chosen, and will be illustrated later in Fig. (8). Similarly 
to Fig. (5) singularity behavior in negativity appears in 
the limit R — > and R — > tt in Fig. (6). 

Similarly, in order to further demonstrate and unify 
the dependence of the entanglement on R (0 < R < tt) 
and B, the negativity surface is plotted as a function of 
R and B for a fixed value T = 0.001 in Fig. (7). This 
surface plot shows the field dependence of the interesting 
valley behavior of the negativity, with side plateaus, as a 
function of R. It is clear that at this fixed value of T = 
0.001, the negativity valley (symmetric in R about R = 
tt/2) widens with increasing temperature, with the side 
plateaus essentially remaining at their low B maximum 
of N = 0.5. The independence of the two symmetric 



Finally we consider hyperbolic version of long-range in- 
teraction in SCM model, which is defined with exchange 



interaction J(R) — Jo/sinh R [46H52| . Here we set 
Jo = 1 and 9 = R for simplicity. The eigenvalues for 
this type of SCM interaction are given by 



(14) 



Using the eigenvalues in Eq. (14) and the eigenvectors 
in Eq. (2) (with Eq. (3)) the density matrix p can be 
constructed. For the SCM interaction type of J(R) = 
1 / sinh 2 R, the partial transpose of p shown in Eq. (4) 
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has matrix elements given in Eq. (15) as follows: 





1 3B 

an = — 

Zj 


(15a) 


a-21 


= >- eT > 


(15b) 


a.33 


= 4(l + e*) 


(15c) 


055 


-> + .*> 


(15d) 


^65 


-id-*) 


(15e) 




1 3B 

&66 = 


(15f) 



where rj = exp {—{B - V2csch[i?] 4 ^sinh [R] 4 ) /2T) , £ = 
exp(-(5 + % /2/^sinh[i?] 4 )/2T). The partition function 
for this case is simply given as 

Z = 2cosh[g] + (l + e B / T )(e + r ? ) . (16) 

Similarly, the negativity follows from summation of the 
negative eigenvalues of the partial transpose density ma- 
trix Eq. (4) with elements as in Eq. (15). 

The negativity for this J(R) = 1/ sinh 2 R case is plot- 
ted versus R (0 < R < it) for different temperature T 
values at fixed magnetic field value B in Fig. (9a), and 
for different magnetic field values at fixed temperature T 
in Fig. (9b). It can be seen from Fig. (9a) that, at a fixed 
value of B — 1, the negativity attains its maximum value 
in a plateau for low temperature, dropping suddenly to 
zero with increasing R at a critical R c value, but it also 
drops to zero in the limit R—>0. In this R — > limit, the 
negativity has a singularity. The value of negativity on 
the plateau and the plateau width in R both vary with 
T, with the lowest value of T = 0.001 exhibiting satura- 
tion at N = 0.5. Fig. (9b) also shows plateau behavior of 
negativity at small R with a sudden drop as R exceeds a 
threshold R c . Similarly to Fig. (9a), in the limit R — > 0, 
a singularity appears in negativity. Negativity also drops 
to zero in the this limit. However, for the fixed value 
T = 1 chosen none of the plateau values for N attain the 
maximum, with the value at the plateau increasing with 
increasing magnetic field B. The plateau width in R is 
essentially independent of B for the range considered. 

In order to further demonstrate and unify the depen- 
dence of entanglement on R and T, the negativity sur- 
face is plotted as a function R (— tt < R < tt) and T for 
a fixed value of B = 1.5 in Fig. (10). This surface plot 
further illustrates the interesting plateau and critical be- 
havior of the negativity. It can be seen from Fig. (10) 
that for increasing \R\ the negativity suddenly drops to 
zero at a critical R c value, but it also drops to zero in 
the limit R — > 0. In this R — > limit, the negativity has 
a singularity since R values are chosen between interval 
—tt < R < it. This Fig. (10) clearly demonstrates that 
the (symmetric positive and negative) critical R c values 
are essentially independent of temperature T, although 
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FIG. 9. Negativity as a function of R for the two-spin fer- 
romagnetic XY model in the case J{R) = 1/ sinh 2 7? (a) for 
selected temperatures at fixed B = 1, (b) for selected mag- 
netic fields at fixed T = 1. 




FIG. 10. Surface plot of the negativity for the two-spin fer- 
romagnetic XY system in the case J(R) = 1/ sinh 2 R as a 
function temperature T and interaction parameter R at fixed 
B = 1.5. 



the actual value of the negativity on the entanglement 
plateau does decrease with increasing T . 

Finally, to further demonstrate and unify the depen- 
dence of entanglement on R and B, the negativity is plot- 
ted as a function R (—it < R < it) and B for a fixed value 
T = 0.1 in Fig. (11). Again plateau and critical behavior 
is apparent. It can be seen from Fig. (11) that the con- 
stant negativity on a plateau suddenly drops to zero at a 
critical R c value. Here the (symmetric positive and neg- 
ative) R c show dependence on B, whereas — aside from 
the small B regime — the value of the negativity on the 
plateau is essentially saturated at N = 0.5. Similarly to 
Fig. (10), in the limit R — > 0, a singularity appears in 
the negativity since R values are chosen between interval 
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FIG. 11. Surface plot of the negativity for the two-spin ferro- 
magnetic XY system in the case J{R) = 1/ sinh 2 R as a func- 
tion magnetic field B and interaction parameter R at fixed 
T = 0.1. 



—it < R < it. In the small B regime the critical be- 
havior disappears and the negativity decreases smoothly 
as the symmetric plateaus disappear for decreasing B or 
increasing R. 

To summarize, in this final subsection, we have inves- 
tigated the behavior of the entanglement of a two-spin 
XY system in the case of the long-range SCM interac- 
tion J(R) = l/sinh 2 i?. We have also discussed the de- 
pendence of the negativity on interaction parameter R, 
magnetic field B and temperature T. We have demon- 
strated the symmetric behavior about R = (with singu- 
lar behavior as R — > 0). Critical plateau behavior is seen 
over a wide parameter regime along with the associated 
entanglement sudden death. 



types of interactions. Our numerical results show that 
in the presence of the SCM type interactions character- 
ized by an interaction parameter R, thermal entangle- 
ment between spin qubits has a rich behavior dependent 
upon R, the temperature T and the applied magnetic 
field strength B. Indeed, we have found that there are 
entanglement plateaus which saturate at the maximum 
negativity N = 0.5 for significant regions of the parame- 
ter ranges of R, T and B. Aside from specific examples, 
in general the entanglement plateaus are characterized by 
a critical distance R c and moving off a plateau through 
R c demonstrates a critical change and entanglement sud- 
den death. Clearly the effect of SCM type interactions 
on the resource of entanglement provides a rich source 
of behavior, with maximum entanglement existing over 
significant parameter regimes and critical switching. 

The plateau behaviors in entanglement formation of 
two-spin given by Eq. (1) for all types of SCM interac- 
tions are very interesting. The physical mechanism of 
entanglement plateaus probably may be long-range in- 
teractions. Indeed, while the plateau behavior in en- 
tanglement does not appear for short-range interactions, 
here, we can see that all long-range interactions with an 
inverse-square and its trigonometric and hyperbolic vari- 
ants lead to different plateau behavior in the entangle- 
ment. Therefore, we can suggest that the plateau be- 
havior of entanglement in these model emerge from the 
long-range character of the interactions. In all case for 
the long-range SCM interactions, entanglement between 
spins still subsists throughout plateau up to critical dis- 
tance although spins are spatially separated from each 
other unlike the system with short-range interactions (67r - 
[7fj| . We believe that these interesting results can be valu- 
able for researcher in this area of physics. 
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